Theory of Extrasolar Giant Planet Transits 
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ABSTRACT 

We present a synthesis of physical effects influencing the observed lightcurve of an extrasolar 
giant planet (EGP) transiting its host star. The synthesis includes a treatment of Rayleigh scat- 
tering, cloud scattering, refraction, and molecular absorption of starlight in the EGP atmosphere. 
Of these effects, molecular absorption dominates in determining the transit-derived radius R, for 
planetary orbital radii less than a few AU. Using a generic model for the atmosphere of EGP 
HD209458b, we perform a fit to the best available transit lightcurve data, and infer that this 
planet has a radius at a pressure of 1 bar, R\, equal to 94430 km, with an uncertainty of ~ 500 
km arising from plausible uncertainties in the atmospheric temperature profile. We predict that 
R will be a function of wavelength of observation, with a robust prediction of variations of at 
least ±1% at infrared wavelengths where H 2 opacity in the high EGP atmosphere dominates. 

Subject headings: planetary systems — stars: individual (HD209458) 



1. Introduction 

Measurements of the diminution of starlight 
during transit of a planet across the disk of a star 
provide an almost direct means to detect extraso- 
lar giant planets (EGPs) with orbital inclinations 
close to 90°. When coupled with measurements 
of the radial velocity variation of the orbited star 
during motion about the common barycenter, the 
mass M of the planet can be measured, and the 
radius R of the planet can be deduced from the 
depth of the transit lightcurve. 

To date, only one transiting extrasolar planet 
has been observed: HD209458b (Charbonneau et 
al. 2000; Henry et al. 2000). A high-quality com- 
posite transit lightcurve has been obtained using 
the STIS spectrograph on the Hubble Space Tele- 
scope (HST; Brown et al. 2000), and a model fit 
to the lightcurve and radial velocity data yields 



the following results: inclination i = 86.68 ± 
0.14°, mass A/=0.69 Mj (where Mj is the mass of 
Jupiter), and radius R = 1.347 ± 0.060i?j (where 
Rj is the radius of Jupiter). Thus, HD209458b has 
been confirmed as a genuine hydrogen-rich EGP 
(Burrows et al. 2000). Spectroscopic radial veloc- 
ity data for HD209458 give a precise value for the 
planet's orbital period, P=3. 524738 days. 

Brown et al. (2000) modeled planet HD209458b 
as a uniform occulting disk of radius R. However, 
as Seager and Sasselov (2000) first pointed out, 
the value of R for a real planet will be a func- 
tion of wavelength, depending on the transmissive 
properties of the planet's atmosphere, as well as 
on other properties of the atmosphere, such as the 
location of dense cloud layers. In this paper, an 
unsubscripted R will denote such a wavelength- 
dependent radius, while subscripted Rs will de- 
note values of the radius at a fixed level in the 
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planet's atmosphere. Specifically, for purposes of 
comparing the inferred values of R and M with 
theoretical models of EGPs of specified age t, one 
should relate the inferred R to the radius of the 
planet at a specific fiducial pressure, as is done 
for Jupiter, where Rj is customarily expressed as 
Rj = 71492 ± 4 km, the equatorial radius at a 
pressure of 1 bar (Lindal et al. 1981). The pur- 
pose of the present paper, then, is to fit the HST 
lightcurve of Brown et al. (2000) to explicit atmo- 
spheric models of HD209458b, in order to derive 
the EGP's radius at 1 bar pressure, which we will 
call R\. Along the way, we obtain further predic- 
tions of the variation of R with wavelength, over 
a broader range of wavelengths than in the analy- 
sis of Seager and Sasselov (2000). Our model for 
the atmospheric structure is a generic one, but it 
differs in some respects from that of Seager and 
Sasselov (2000). 

The transit lightcurve depends on (a) Raylcigh 
scattering of light from the host star, (b) refraction 
of the stellar surface brightness distribution, and 
(c) the slant optical depth r through the planet's 
atmosphere, as determined by molecular opacity 
and clouds. All of these effects depend in turn 
upon the atmospheric pressure (P) vs. tempera- 
ture (T) profile, and upon the surface gravity g. 

In the following, we consider the P vs. T profile, 
effects (a) through (c) , and then present results for 
the relation of R as a function of wavelength A and 
the best-fit results for i?i for HD209458b. 

2. Atmospheric Model 
2.1. P-T Profile 

Our philosophy in constructing the P—T profile 
of our baseline model for HD209458b is that this 
profile should be representative for the planet's 
albedo class (either Class IV or V), as defined in 
Sudarsky, Burrows, and Pinto (2000, SBP). The 
specifics of this profile are not important as long 
as the basic molecular composition of the atmo- 
sphere is respected and the mapping between pres- 
sure and areal mass density is correct for a given 
gravity. The surface gravity of HD209458b is mea- 
sured to be close to the Earth's value. Hence, 
we used the [T off =1270 K/gravity=10 3 cm s~ 2 ] 
model for Class IV EGPs, similar to that found in 
SBP. (Note that the actual P-T profile will be 
a function of dynamic processes that redistribute 



heat from pole to equator and to the night side, 
processes that are difficult to model [Guillot and 
Showman 2001].) 

Figure 1 shows the nominal P — T profile used 
in the model (solid heavy line on left side), along 
with high-temperature and low-temperature ver- 
sions used to explore the sensitivity of the results 
to the profile. At P = 1 mbar, the difference be- 
tween the high-temperature and low-temperature 
profiles is about 500 K, or about 50 percent of 
the nominal temperature. As we discuss below, 
the transit radius of HD209458b is primarily de- 
termined by opacity sources in the vicinity of P 
~ 10 mbar. In §7 we discuss the sensitivity of 
the transit radius to the alternative temperature 
profiles. 

The right-hand side of Fig. 1 shows adiabatic 
P — T profiles for various interior models of the 
planet, to be discussed in §8. 

2.2. Clouds and Condensates 

We predict the altitude at which clouds of var- 
ious condensable species form using a code sum- 
marized in Burrows et al. (2001) and Marley et al. 
(1999). Vapor pressure relations for the rocky con- 
densates are from Lunine et al. (1989). Above the 
altitude at which a given condensate first appears, 
growth rates for particles and droplets are calcu- 
lated using analytic expressions (Rossow 1978). 
The cloud model assumes that the atmospheric 
thermal balance at each level is dominated by a 
modal particle size that is the maximum attain- 
able when growth rates are exceeded by the sed- 
imentation, or rainout rate, of the particles. In 
convective regions, equating the upwelling (con- 
vective) velocity and the sedimentation velocity 
sets the particle size. The amount of condensate 
at each altitude within the cloud is given by the 
vapor pressure at that level multiplied by a "super- 
saturation factor" usually set to 0.01 by analogy 
with terrestrial clouds. The particle size and the 
mass density of the condensate at each level thus 
determine the number density of particles. 

The enstatite cloud layer is potentially the most 
important for affecting the value of R. Scatter- 
ing and extinction cross sections for this major 
cloud forming species were obtained for the com- 
puted particle sizes by a full Mie theory (SBP). 
The optical properties of enstatite were taken from 
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Dorschner et al. (1995). 

Ackcrman and Mar ley (2001) develop a model 
of cloud formation that extends the foregoing 
to include a more realistic rainout prescription. 
However, as shown in Fig. 2, the major cloud 
forming species for HD209458b, enstatite, does 
not contribute opacity in the right altitude range 
to have a significant effect on the transit pro- 
files. More refractory cloud- forming species, e.g., 
aluminum silicates, occur deeper in the atmo- 
sphere. Less refractory cloud formers such as wa- 
ter and sulfur-bearing species would condense out 
at higher altitude, but the atmosphere is too warm 
for these species to occur as clouds. Therefore, for 
HD209458b, our cloud model is more than ade- 
quate. It is possible that very minor cloud form- 
ing species that are slightly less refractory than 
enstatite would put some cloud opacity at mod- 
estly higher altitudes, but their smaller abundance 
relative to enstatite would proportionately dimin- 
ish their effect. We therefore argue that, for this 
particular extrasolar planet and those with similar 
effective temperatures, cloud opacity is not signif- 
icant in determining the transit radius. Adopt- 
ing the Ackerman and Marley (2001) prescription 
would likely result in an even less significant role 
for the enstatite clouds, since such a model would 
result in rainout of more condensate and lead to a 
less optically-thick cloud. 

3. Rayleigh Scattering 

Rayleigh scattering was treated in an approxi- 
mate manner. As we demonstrate below, Rayleigh 
scattering has only a minor effect on the value of 
R, so an approximate treatment is sufficient. The 
effect of Rayleigh scattering (or any other scat- 
tering) is complex, because stellar photons inci- 
dent on the planet's atmosphere in a given pen- 
cil with incident intensity / are partially removed 
from the pencil and scattered into different solid 
angles. Thus, if only conservative Rayleigh scat- 
tering were occurring, the limb at R would be 
defined by the radius at which a high probabil- 
ity of such removal occurs. However, photons are 
scattered into the beam to the observer as well 
as removed. Rather than treat the full three- 
dimensional problem of Rayleigh scattering in a 
spherically-stratified atmosphere, we replaced the 
atmosphere with a series of slabs located in a plane 



containing the center of the planet and orthog- 
onal to the star-observer line. In the following, 
we will denote the two-dimensional vector separa- 
tion of a point in this plane from the projected 
planetary center by r, and the scalar value by r. 
The three-dimensional vector position of an at- 
mospheric point from the planetary center will be 
denoted by r". 

Each slab, located at a two-dimensional radius 
r from the projected planet center, has a Rayleigh- 
scattering optical depth tr given by 
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where N is the number density of molecules in 
the atmosphere, and ctr is the Rayleigh-scattering 
cross-section per molecule, given by 
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where v is the refractivity (refractive index minus 
1) of the gas (Chandrasekhar 1960). Since v <C 1 
and v cx N, or is a function only of A and the gas 
composition. 

We wrote a Monte Carlo scattering code to in- 
vestigate the effects of Rayleigh scattering in the 
planet's atmosphere. The code follows every pho- 
ton as it travels through a plane-parallel slab of a 
given optical thickness. Incident photons can ar- 
rive at the top of the slab from any direction on an 
imaginary hemisphere. Inside the slab, after every 
scattering event, a new photon direction in three 
dimensions is calculated, based on the Rayleigh- 
scattering phase function. Photons are followed 
until they emerge from either side of the slab, and 
the total path traveled, in units of optical depth, 
is tabulated. When the photon finally emerges 
from the bottom of the slab it is placed into a bin 
that corresponds to its final direction. Upon com- 
pletion, the total number of photons in each bin 
is obtained, along with the average path traveled 
per photon in that bin. 

Our code was tested against the analytical re- 
sults of van de Hulst (1974) and Chandrasekhar 
(1960). To within the noise in the Monte Carlo 
simulations (~ 1%), we were able to match van de 
Hulst 's results for an isotropic distribution of pho- 
tons and Chandrasekhar's for isotropic and "pencil 
beam" distributions at a variety of incident angles. 
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We ran our simulations for radiation at normal 
incidence to the slab for a variety of different op- 
tical thicknesses, logarithmically spaced from 0.01 
to 28. If we assume an imaginary observer looks at 
normal incidence from the other side of the slab, 
the observer will see photons that pass through the 
slab unobstructed, as well as those that are multi- 
ply scattered back into the beam and emerge nor- 
mal to the surface. Of most interest to us in this 
situation are those photons that, although scat- 
tered, emerge normal to the slab surface and ulti- 
mately reach the observer. These photons create 
a Rayleigh-scattering "glow" from the slab. We 
found that the glow intensity in the direction of 
the observer increases until tr ~ 3.4, as fewer and 
fewer photons are able to pass through the slab 
unscattered. At tr > 3.4, up to tr = 28, the 
greatest thickness we ran, Rayleigh-scattered in- 
tensity decreases more or less exponentially, 
larger fraction of the photons are scattered back 
out the top of the slab, rather than scattering all 
the way through. We fitted an equation to both 
sides of this curve so that the glow intensity could 
be interpolated, and extrapolated to higher tr if 
necessary. We also fitted an equation to the aver- 
age optical path length traveled for these photons, 
using this result to estimate the total molecular- 
absorption optical depth, tm, and cloud optical 
depth, to The total optical depth for photons 
initially incident with impact parameter r is then 
given by t = ir + tm+tc- In practice, for all cases 
that we have investigated, tm becomes large long 
before tr or tc does. Likewise, refractive effects 
are still negligible when tm becomes significant, as 
discussed below. Thus, for purposes of computing 
tm, the photon path can be taken to be a straight 
line through the planetary atmosphere, such that 

/>oo 

T M = 2 / r"dr"a M N(r")/Vr" 2 -r 2 , (3) 

J r 

where cm is the average absorption cross-section 
per molecule for a solar-composition planetary at- 
mosphere. We estimate tc with a similar formula, 
although strictly speaking, it becomes appreciable 
at levels where refraction is important. 

Since in our transit calculations the planet is a 
disk passing in front of its star, we assume the at- 
mosphere of the planet is a flat slab with an optical 
thickness that decreases with increasing distance 
from the planet center. In this way, our results for 



different slab thicknesses can be used directly in 
our planetary atmosphere calculations. 

In our transit simulations, the atmosphere of 
the planet was broken up into many small an- 
nuli, each with its own optical thickness, which 
is dependent on wavelength, as calculated using 
our model for the structure of the planet's atmo- 
sphere. At a given wavelength of light, for every 
optical thickness in the atmosphere, the Raylcigh 
glow intensity was calculated. This intensity re- 
sults in a small additive component to the star 
plus planet signal outside of transit, but in gen- 
eral has no detectable effect on the observed light 
curve. In addition, since the total average path 
traveled for a scattered photon is also known, the 
effects of absorption due to molecular opacity can 
be calculated. At high optical thicknesses, most 
scattered photons that would emerge from the slab 
are actually absorbed by opacity sources in the at- 
mosphere. 

4. Refraction 

The theory used to compute refractive effects is 
essentially identical to that of the standard theory 
for occultations of stars by planetary atmospheres 
(Hubbard, Yelle, and Lunine 1990). In the follow- 
ing, we let I(r) denote the photon intensity within 
a differential solid angle whose coordinates are la- 
beled by the two-dimensional vector r measured in 
the plane of the sky, from the center of the planet. 
The intensity I(r) of the observed image of star 
plus planet is then given by the mapping 

I{v)=I'{v')e-^\ (4) 

where I'(r') is the stellar surface brightness dis- 
tribution, t is the total optical depth integrated 
along a ray path with impact parameter r, and 
r' and r are two-dimensional vectors in a plane 
normal to the propagation direction, marking the 
starting point of a ray on the stellar surface and 
its closest- approach position in the planet's atmo- 
sphere, respectively. 

The mapping from r' to r is obtained by com- 
puting the total phase shift $ imposed by the rc- 
fractivity distribution on a photon with impact pa- 
rameter r, 



= (4tt/A) / r"dr"v(r")/\/r" 2 -r 2 , (5) 



4 



then computing the two-dimensional bending an- 
gle a. according to ol = (A/27r)V<I>, where the gra- 
dient is taken in the two-dimensional plane. 

Figure 2 shows the calculated values of tm, tr, 
and cloud optical depth, tq , for our best-fit plane- 
tary model, as discussed in §6. Figure 2 also shows 
the atmospheric radius where refraction begins to 
be important. Specifically, the shaded region in 
Fig. 2 shows the difference between r and r' at a 
level where the difference, 

r — r' = Da, (6) 

becomes equal to 500 km, or about one atmo- 
spheric scale height (D is the distance from the 
star to the planet). The upper edge of the shaded 
region corresponds to r and the lower edge to r'. 
Since a varies exponentially with — r, consider- 
able ray bending occurs for impact parameters be- 
low the shaded region. However, tm is always so 
large in this region that refraction is unimportant 
for defining R, for values of D similar to that of 
HD209458b. 

We do not include general-relativistic ray bend- 
ing in the theory presented in this paper, as it 
is unimportant for the parameters of HD209458b. 
However, it is straightforward to include gravita- 
tional lensing. One simply adds to the refractive 
a discussed above, an additional term for general- 
relativistic ray bending, «gr = AGM/rc 2 (where 
G is the gravitational constant, M is the planet's 
mass, and c is the speed of light). The total bend- 
ing angle used in Eq. (6) is then a + acR- 

5. Gaseous Opacities 

The primary gaseous absorptive opacity sources 
in the atmospheres of hot EGPs include H 2 , H 2 0, 
CH4, CO, and the important alkali metals, Na 
and K. We take the temperature- and pressure- 
dependent opacities from theoretical and experi- 
mental data referenced in Burrows et al. (f997) 
and Burrows et al. (2001). In the near infrared 
(~1 to 2.6 //m), absorption by H 2 molecules fig- 
ures most prominently, with strong ro-vibrational 
bands centered at ~0.95, 1.15, 1.4, 1.85, and 2.6 
fim. The strong pressure-broadened resonance 
lines of neutral Na and K - the strengths of which 
depend on the level of ionization by stellar irra- 
diation - appear prominently at ^0.59 and 0.77 
fim, respectively. The dominant carbon-bearing 



molecule is a function of both temperature and 
pressure. At lower temperatures, CH 4 is domi- 
nant, but CO overtakes CH4 at higher temper- 
atures (~950 K at 0.1 bar, -1100 K at 1 bar). 
Hence, the strengths of CH4 features (1.4, 1.7, 
2.2 /im) and CO features (1.2, 1.6, 2.3 fim) arc 
highly temperature-dependent. Finally, an im- 
portant continuous opacity source in cloud-free 
EGPs is H 2 -H 2 collision induced absorption at 
high pressures and temperatures (Zheng and Bo- 
rysow 1995). 

All of the opacity sources mentioned above were 
used to calculate cm, and then incorporated in Eq. 
(3); the results are plotted in Fig. 2. 

6. Transit Lightcurve 

The next step in calculating a transit lightcurve 
was to synthesize images of the two-dimensional 
distribution of starlight around the planet. We 
created a synthetic square aperture of size 361 x 
361 pixels (1 pixel = spatial scale of 700 km), cen- 
tered on the planet. The star was taken to be a 
disk with a pixel intensity of 1 at its center, with 
a prescribed darkening law to the limb. At each 
pixel in the array, the pixel intensity was calcu- 
lated using the calculated values of the various rs, 
adding the Rayleigh-scattered component into the 
beam and the refracted stellar component, from 
each pixel on the stellar disk (including contribu- 
tions from virtual pixels on the part of the stel- 
lar profile outside the synthetic aperture). The 
total pixel sum over the aperture was then calcu- 
lated, with and without the planet present, allow- 
ing us to calculate the total intensity subtracted 
and added by the planet. 

We then synthesized a transit lightcurve, incor- 
porating all the physical effects described above, 
and using most of the fitted parameters of Brown 
et al. (2000): i, stellar radius, and orbital radius 
and period. The lightcurve was obtained by aver- 
aging over the bandpass of the Brown et al. (2000) 
experiment, weighted by a blackbody distribution 
for the effective temperature of HD209458. The 
value of R\ for HD209458b was adjusted in our 
model until we matched the depth of the theoret- 
ical transit lightcurve to the depth of the compos- 
ite observed lightcurve, as shown in Fig. 3. The 
match of our model, which has only the adjustable 
parameter Ri, to the data is quite good. We ob- 
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tain #1 = 94430 km for HD209458b, for the nom- 
inal P — T profile shown in Fig. 1. For the cold 
profile, we obtain R\ = 95000 km, and for the hot 
profile, Ri = 94260 km. 

The flux from the star was computed us- 
ing limb-darkening coefficients from Van Hamme 
(1993). The two equations used were the nonlin- 
ear logarithmic and square-root laws. These laws 
are (logarithmic): 

= l-A(l- l £)-Bn\n{ l j), (7) 

and (square root): 

= 1 - C(l - /i) - £>(1 - y/P), (8) 

where 1(1) is the specific intensity at the center of 
the stellar disk, A, B, C, and D are wavelength- 
dependent constants, and fi is the cosine of the 
angle between the line of sight and the emergent 
intensity. 

The curve shown in Fig. 3 is for the logarithmic 
law; we found that use of the square root law in- 
stead changed the lightcurve by less than 1 part in 
400. The coefficients used were interpolated from 
the Van Hamme monochromatic coefficient tables, 
for logg = 4.3, T g of 6000 K, and solar metallic- 
ity The very slight mismatch with the data is per- 
haps due to a difference in the calculated stellar 
intensities, which are based on the stellar models 
of Kurucz (1991), and the actual change in stellar 
intensity across the disk. 

Figures 4 and 5 show examples of synthetic im- 
ages and provide more details about the contribu- 
tions of Rayleigh scattering and refraction to the 
shape of a transit lightcurve, which are generally 
negligible for the specific case of HD209458b. In 
Fig. 4a, we show that even at a short wavelength 
of 0.45 /im, the forward-scattered Rayleigh glow 
is almost invisible. In Fig. 4b, we reduce tm by 
a factor 10 3 and stretch the image to make the 
narrow Rayleigh glow annulus more apparent. 

In Fig. 5, we retain the atmospheric structure 
for the nominal model corresponding to Fig. 4(a), 
but we increase the distance D to 10 4 times the 
value for HD209458b, thus bringing the refracting 
layers of the atmosphere into play. According to 
Eq. (6), for r — r' fixed at, say, one scale height 
(<~ 500 km), then if D increases by 10 4 , a must 
decrease by the same factor. Since a is roughly 
proportional to pressure in the atmospheres con- 



sidered here, a decrease in a by this factor corre- 
sponds to displacing the layer responsible for the 
prescribed refraction upward by about nine scale 
heights, to a region where the slant optical depth 
r is negligible. The atmospheric depth where the 
limb of bright refracted starlight terminates is de- 
termined by the abrupt increase of r to apprecia- 
ble values. 

As D increases further, our theory smoothly 
transforms to conventional ray-optical stellar oc- 
cultation theory in the limit r = (Hubbard, 
Yclle, and Lunine 1990), appropriate to observa- 
tion of the passage of a planet in front of a star of 
smaller apparent angular size than the planet. 

Except for such conventional stellar occulta- 
tions, we are unlikely to observe a transit domi- 
nated by refractive effects, as the probability of a 
transit by an EGP at an orbital radius of hundreds 
of AU is exceedingly small. 

7. Variation of R with A 

As discussed above, the relation R(X) is al- 
most entirely determined by molecular absorption 
features. Seager and Sasselov (2000) discussed 
a possible dramatic variation of R(X) due to al- 
kali absorption features in the visual wavelength 
bands. We predict dramatic effects at infrared 
wavelengths as well, due to strong features of H 2 0; 
however, detection of infrared variations in R(X) 
may require observations above the Earth's atmo- 
sphere. 

Figure 6 shows, for HD209458b parameters, 
the i?(A) relation predicted by our three P — T 
models. This relation is computed by evaluating 
binned averages of opacities (for specified temper- 
atures) at 500 individual frequencies spaced across 
several xlO 4 original frequencies. Note that the 
variation of R vs. A increases with T at some 
wavelengths, and decreases at others. Note also 
that though a flattened temperature gradient in 
the upper atmosphere will smooth out the reflec- 
tion/emission spectrum of the planet itself, the 
variation i?(A) with wavelength will not be sim- 
ilarly flattened. For example, even though the 
Na/K alkali metal features in the planet's spec- 
trum may be smoothed by stellar irradiation, the 
transit size will still vary appreciably across these 
features as long as Na/K are not ionized. 

Because the variations with wavelength are 
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quite rapid in some wavelength intervals, it may 
be possible to strategically choose wavelengths of 
observations which will span these rapid varia- 
tions and still be close enough together to allow 
the limb-darkening of HD209458 to be removed by 
fitting a smooth model. Note that R is typically 
about 2000 km larger than R\ , and thus pertains 
to atmospheric layers which are at pressures near 
<~ 10 mbar. 

8. Conclusions 

Referring back to Fig. 1, note that an adia- 
bat in a solar-composition object with a mass of 
0.69 Mj, which would yield a value of R\ compat- 
ible with the one determined here, must have a 
significantly lower specific entropy than the atmo- 
spheric layers heated by the star. It follows that 
there must be a significant region in the planet, 
possibly spanning pressures from <~ 10 to <~ 10 4 
bar, where the T vs. P relation must be substan- 
tially subadiabatic or even isothermal. Detailed 
thermal evolution models for this layer remain to 
be calculated. 

We have shown that apparent variations of 
at least ±1% in the radius of giant planets oc- 
cur as a function of wavelength. This variation 
is potentially discernable with the next genera- 
tion of transit-observing spacecraft in Earth or- 
bit, provided that they possess the capability of 
multi-wavelength observation with sufficient spec- 
tral resolution. We therefore recommend that 
proposers or designers of such missions consider 
the possibility of transit-based probing of extra- 
solar planet atmospheric composition, via multi- 
wavelength observations. 

Finally, as is suggested by Fig. 6, there could 
be considerable variation in i?(A) at ultraviolet 
wavelengths. Indeed, as is well known (Smith and 
Huntcn 1990), the location of the level where r <~ 1 
in our own giant planets is a strong function of ul- 
traviolet wavelength, and in the UV is typically 
at much lower pressures than 1 /ibar. Observa- 
tion of an EGP transit at UV wavelengths is, in 
many experimental aspects, equivalent to a solar- 
occultation UV experiment in our solar system, 
and may have similar diagnostic power for chemi- 
cal composition of the high planetary atmosphere. 
To be sure, Jupiter has a rather warm strato- 
sphcre/mesosphere, with photochemical aerosols, 



so the UV radius could be a sensitive function of 
the external boundary conditions within which an 
EGP planet finds itself. 
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Fig. 1. Heavy line on the left shows the nominal 
atmospheric pressure-temperature profile used in 
these calculations; light solid line shows the low- 
temperature version and open-dotted line shows 
the high-temperature version. Dashed lines show 
interior H-He adiabats, with the corresponding 
entropy (in units of Boltzmann's constant per 
baryon) and total radius R\ for a planet mass of 
0.69 Mj. Heavy dashed line shows the approx- 
imate interior adiabat that matches the transit- 
derived value of i?i . 
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Fig. 2. — Radius vs. optical depth for tm (dashed), 
tr (dotted), and r c (dot-dashed), for the nominal 
P-T profile and for parameters of HD209458b. 
Below the shaded zone near 1 bar pressure, refrac- 
tion is also significant. Molecular opacity profiles 
are plotted for wavelengths spanned by the obser- 
vations of Brown et al. (2000). 
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Fig. 3. — Synthetic transit lightcurve for our 
model (solid curve), compared with composite 
lightcurve of Brown ct al. (2000). The three P — T 
profiles give indistinguishable solid curves. 




Fig. 4. — (a) Synthetic image of a transit of 
HD209458b (pixel scale is from to 1, with 1 corre- 
sponding to the intensity at the center of the stel- 
lar disk), with the planet at an orbital position 5° 
from inferior conjunction, at a wavelength of 0.45 
[Mm. Close scrutiny will show a faint Raylcigh- 
scattering ring around the planet's limb exterior 
to the stellar limb, (b) An image of the same ge- 
ometry, but with molecular opacity reduced by a 
factor 1000, and with the pixel scale stretched over 
the range to 0.1. The Rayleigh-scattcring ring is 
thus wider and brighter. 




Fig. 5. — Image of the same atmosphere as in Fig. 
4(a) (pixel scale is again from to 1), but with the 
distance D from the planet to the star increased 
from the value for HD209458b to a value 10 4 times 
greater. The bright limb outlining the planet's 
disk is far brighter than the Rayleigh "glow" , and 
is produced by refractive mapping of the stellar 
disk's intensity distribution. 
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Fig. 6. — Predicted variation of transit radius R 
with wavelength (heavy curve, nominal P — T pro- 
file; light curve, cold P — T profile; open dots, hot 
P — T profile; see Fig. 1). The right-hand scales 
show, in percent, the variation of R with respect to 
R\ and with respect to Rt, an ^"average" transit 
radius in the visual wavelength band, adopted as 
Rt = 96500 km = 1.35i?j. At wavelengths where 
slant optical depth is high, R is larger. "Absorp- 
tion" features thus appear upside-down on this 
plot. Prominent features are labeled with the re- 
sponsible molecule. 
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